A computational investigation towards substitution effects on 8π electrocyclisation of conjugated 1,3,5,7-octatetraenes

A computational investigation using M06-2X/6-31+G(d) method is reported for the substitution effects on 8π electrocyclisation of conjugated octatetraene. This systematic study describes the mono- and di-substitution effect across the 1,3,5,7-octatetraene skeleton. A general preference of the outward substitution over the inward, at C1 position of the monosubstituted system is observed. However, mesomerically electron donating group (–NH2 and –OH) display an opposite effect with respect to secondary orbital interaction (SOI) between the lone pair on the substituent and the orbital. A comparative evaluation on the computed activation energies for the 1-, 2-, 3-, and 4-monosubstituted system showed an insignificant impact on the rate of the reaction, in contrast to the electrocyclic ring closure of the unsubstituted compound. Computations of disubstituted system are more pronounced, where a remarkable acceleration is observed for 2-NO2–7-NO2 substituted octatetraene at 4.9 kcal mol−1, and a noticeable deceleration for 4-CH3–5-CH3 substituted octatetraene at 25.4 kcal mol−1 from the parent molecule, 17.0 kcal mol−1. A visible accelerated effects are commonly exhibited by the substitution on the terminal double bonds (C1, C2, C7, and C8), that are 1,2-, 1,7-, 1,8-, and 2,7-patterns, in regard to the greater orbital interaction for the new σ-bond formation. Despite the unfavourable steric clashes of the substituents in the 1,8-system, an apparent reduction in the energy barrier up to 7.4 kcal mol−1 is computed for 1-NH2–8-NO2 system from 17.0 kcal mol−1. This is due to the synergistic effect of the electron donor and electron acceptor, enhancing the stability of the transition structure. The electrocyclic ring closure involving vicinal substitution patterns, such as 1,2-, 2,3-, 3,4-, and 4,5-systems are critically dominated by steric crowding between the adjacent functional groups. In certain cases of the 1,2-substituted system, a noticeable accelerated effects are found for 1-NH2–2-NH2-substituted compound (9.7 kcal mol−1) due to an increased in electronic density on the substituted terminal double bond (C1–C2), hence favouring the formation of the new σ-bond.


Introduction
Carbocycles 1 are structurally diversied motifs, which are ubiquitously present in bioactive natural products and cytotoxic molecules.
In general, electrocyclisation 11 is a powerful methodology to construct polycyclic scaffolds with highly efficient atom economy. 12It proceeds with an excellent and predictable regioand stereocontrol, governed by Woodward-Hoffmann's rules. 13ccording to these rules, electrocyclic ring closure of conjugated octatetraenes with 8p system is a thermally induced conrotatory process, which involves spontaneous rotation of the terminal olen bonds of a tetraene.Particularly, the steric effect, electronic properties and torquoselectivity 14 (inward/outward rotation) of the substituents located at the C-atoms subjected to electrocyclic transformation (ring-opening or closure) play a vital role in orchestrating the observed stereoselectivities.For the conrotatory electrocyclisation of 1-substituted 1,3,5,7-octatetraenes, a pronounced steric effect of the substituents inuences the stereochemical outcomes as reported by Houk and colleagues using ab initio molecular orbital theory.14b In addition, electronic properties of the substituents renders an insignicant effect on stereoselectivities for the octatetraene cyclisation due to dearth of rotational preference caused by the inward and outward substituents located at the terminal Catoms involved in the bond formation.These ndings are in sharp contrast to the conrotatory ring-opening of 3-substituted cyclobutenes, in which, the stereoselectivities are primarily mediated by the electronic nature of the substituents.
Building on the pioneer work championed by Houk et al., 14 Cossío and co-workers exploited the application of density functional theory (DFT) to further elaborate the electrocyclic ring closure for variously 1-substituted octatetraenes (Fig. 1C(I)). 15A detailed investigation suggests complete periselective electrocyclic ring closure of octatetraenes via Möbius aromatic transition structure.With a couple of exceptions (-NH 2 and -NH 3 + substituents), lower barriers have been found for the transition structures which have substituents at the outward positions and this is independent of the electronic nature of the substituents.The lowest activation energy (11.8 kcal mol −1 ) is computed for inward NH 3 + -substituted octatetraene, which is considerably lower than the measured activation energy for the unsubstituted octatetraene i.e., 17.0 kcal mol −1 . 16In addition, rate of reactions for 1,8-dimethyl substituted octatetraene isomers have also been computed, which require substantially higher barrier compared to that of the parent compound.Later, the groups of Schreiner and Suffert applied a computational approach to explore the origin of high torquoselectivities and complete diastereoselectivities observed during the competing synthesis of cyclooctatrienes vs. fenestradienes via 8p electrocyclic transformations. 17As expected, Möbius aromatic transition structures having different P-and M-helical topologies led to origin of observed selectivities in the cyclooctatrienes, formed as a result of 8p conrotatory electrocyclisation of tetraenes.However, the corresponding fenestradienes are the products of the 8p-6p electrocyclic cascade reactions.
Owing to synthetic utilities and increasing research interests towards 8p-6p electrocyclisation cascade, an inclusive and efficient computational investigation dealing with terminally CH 3 -substituted 1,3,5,7-octatetraenes has been elegantly reported by Patel and Houk (Fig. 1C(II)). 18The inuence of steric crowding at terminal positions of the tetraene is directly linked with the reaction thermodynamics.Formation of corresponding 1,3,5-cyclooctatrienes via 8p electrocyclisation from the un-and monosubstituted tetraenes proceeds rapidly however, 6p electrocyclization leading to highly diastereoselective bicyclic product is relatively slow and reversible.In addition, initial 8p electrocyclisation for di-and trisubstituted reactants is comparatively less exergonic due to destabilization rendered by the steric effect.This leads to an endergonic 8p electrocyclic product for tetrasubstituted tetraene however, the overall reaction involving this cascade still remains exergonic and yields a bicyclic product containing two vicinal quaternary centers.
Here, a systematic and succinct DFT study is reported which delineates the mono-and disubstitution effect across the skeleton of 1,3,5,7-octatetraenes (Fig. 1D).Such an effect has already been reported for thermal conrotatory 4p electrocyclic ringopening of substituted cyclobutenes 19 and thermal disrotatory 6p electrocyclic ring-closure of 1,3,5-hexatrienes. 20o start with, a selection of commonly used substituents ranging from electron donating to withdrawing substituents have been subjected to this study.A remarkable accelerated substitution effect has been noticed that lowers the activation free energy for the 8p electrocyclisation up to 4.9 kcal mol −1 compared to 17.0 kcal mol −1 , experimentally measured for the unsubstituted 1,3,5,7-octatetraenes. 16At the same time, certain disubstituted patterns produce negative effect, leading to a noticeable increase in the activation barriers; the maximum computed to be 23.6 kcal mol −1 .

8p electrocyclisation of unsubstituted 1,3,5,7octatetraenes
At the outset, a detailed comparative benchmarking study was envisaged as both B3LYP and M06-2X methods have been reported for the electrocyclisation of octatetraenes (Table 1). 15,18or an extensive and systematic study of a particular transformation, as it is, computational expense plays a crucial role, as well as, evaluation of calculations with the experimental values.Several basis sets were selected for both methods to nd an appropriate computational approach for further investigation.It is evident from this study that B3LYP method overestimates the barrier for tetraene by ∼2 kcal mol −1 compared to M06-2X method.With the same time, a signicant difference in reaction energies have been observed for both methods.Apparently, B3LYP method fails to accurately compute the reaction thermodynamics by showing approximate thermoneutral behaviour of electrocyclisation while M06-2X method, as expected, describes the exergonic nature of the 8p electrocyclisation of unsubstituted 1,3,5,7-octatetraenes.
Compared to experimental result, M06-2X/6-31+G(d) is the most appropriate basis set for the unsubstituted 1,3,5,7-octatetraenes (entry 2, Table 1).It is worth noticing that an absence of diffused function (entry 1, Table 1) and addition of polarization function to the hydrogen atoms (entry 3, Table 1) result in lower activation barriers compared to experimental value.At the same time, basis set containing both diffused and polarizations functions for heavy and hydrogen atoms offer slightly higher barrier (entry 4, Table 1).The reaction preferably proceeds through the helical shaped TS, 21 which is associated with the conrotatory electrocyclisation of tetraene.In addition, disrotatory movement of molecular orbitals would undergo via deformed boat-shaped TS of the 1,3Z,5E,7-octatetraenes and has signicantly higher activation barrier, as reported by Cossío et al. 15 For the optimum method concluded from this benchmarking study [M0-62X/6-31+G(d)], the activation energy difference between the conrotatory and disrotatory TSs is 28.8 kcal mol −1 , which is substantially higher than the measured activation energy i.e., 17.0 kcal mol −1 . 16More or less, Table 1 Comparative benchmarking study using B3LYP and M06-2X methods to compute the activation energies (DG ‡ , kcal mol −1 ) and reaction energies (DG, kcal mol −1 ) for conrotatory 8p electrocyclic ring-closure of 1,3Z,5Z,7-octatetraene and disrotatory 8p electrocyclic ring-closure of 1,3Z,5E,7-octatetraene (in parentheses)

Entry
Basis set B3LYP M06-2X similar difference in the activation barriers is observed for all the basis sets considered for this study.Further comparative analysis between B3LYP and M06-2X methods reveals a good correlation between the activation barriers for the conrotatory and disrotatory electrocyclisation using different basis sets for B3LYP method (Fig. 2a) and the corresponding M06-2X method (Fig. 2b).There is a strong correlation between the activation and reaction energies for the conrotatory electrocyclisation, however poor correlation was found for the disrotatory process.

8p electrocyclisation of monosubstituted 1,3,5,7octatetraenes
The computational results of thermally allowed conrotatory 8p electrocyclic ring closure of terminally monosubstituted system using M06-2X/6-31+G(d) method are shown in Table 2.The rate of these octatetraene cyclisation is determined by the outward or inward arrangement of the substituent on the system.Cossío et al. has identied a general preference of the outward attachment over the inward, as a consequence of lower steric crowding experienced by the substituent in the transition state. 15Such interaction renders slight stabilization to the outward transition structure, which was reected by the lower computed activation barrier regardless of its electron donating or electron accepting character.Our computational ndings have also shown similar trends.A minor preference of 1.5 to 2.6 kcal mol −1 lower than the computed activation energies for the substituent to be in the inside position are recorded.Consequently, the stereochemical outcomes of these electrocyclic reactions are directed by the steric effects.
On the other hand, installation of mesomerically electron donating group, -NH 2 and -OH group, shows an opposite preference of the torquoselectivity in the electrocyclisation process.Presumably, a secondary orbital interaction (SOI) 15 between the lone pair on the functional group and the localized p * C7;C8 orbital, as shown in Fig. 3, provide extra stabilization to the transition structure.This interaction is favoured by the excellent overlap of the atomic orbital of the inward substituent with the atomic orbital on C7 (Fig. 3a).Moreover, it depends on the capacity of the Fig. 2 The correlation between activation free energies for conrotatory (con) and disrotatory (dis) electrocyclisation of unsubstituted 1,3Z,5Z,7-octatetraene and 1,3Z,5E,7-octatetraene respectively, using different basis sets for (a) B3LYP and (b) M0-62X methods.
electron donor character.Introduction of other mesomerically electron donating groups, -CH 3 and -F, at the terminal position with an inward substitution does not offer a stabilizing SOI to the transition structures (Fig. 3b and e), particularly due to the weaker electron donating effect.Despite the presence of lone pairs on -F, the large electronegativity demonstrates an inductive effect that lowers the p-donor capability towards the p * C7;C8 orbital of the terminal double bond (Fig. 3e).
It is noteworthy that the carbonyl (-CHO) and hydroxyl (-OH) functional groups can be oriented in two possible conformations on the tetraene molecule, as was reported in the computational investigation of substituted 4p system paper. 19ence, lowest energy conformation is used to study the reaction energies.In contrast to the unfunctionalised octatetraene (17.0 kcal mol −1 ), electronic properties of the substituent in a 1substituted system hardly affect the rate of the transformation, ranging from 14.9 to 20.7 kcal mol −1 .The mild steric effects do bring about the minute preference for the outward mode of rotation of the terminal substituent.
We have also computed the activation and reaction energies of the transition structure with functionality located at 3 other possible positions (C2-, C3-, and C4-position) to extensively study the mono-substituent effect on the reactivity of electrocyclization of the tetraene molecule (Table 3).Based on the computations, electronic properties of the functionality at C2position are maximized, while the steric effects are minimized.A much larger reduction in the activation energies is observed, contrary to the other monosubstituted systems.At C3and C4-monosubstituted octatetraenes, the energies are practically indistinguishable from the terminally monosubstituted system.The reasons may be attributed to the inconsequential electronic interaction of the substituent due to an indirect interaction with the terminal double bonds that are involved in the formation of the new s-bond.
Particularly, incorporating a strong electron-withdrawing group (-CHO and -NO 2 ) at these positions (C2, C3 or C4) reveals an increase in the rate of the tetraene ring closure.A noticeable reduction is shown especially by the substitution at position C2, where the electronic properties are maximized and steric congestion are minimized.
The reaction of 2-CHO (Table 3, entry 5) and 2-NO 2 (Table 3, entry 7) system are accelerated to 11.4 kcal mol −1 and 10.2 kcal mol −1 , respectively, from 17.0 kcal mol −1 as reported Fig. 3  Table 3 Activation energies (DG ‡ , kcal mol −1 ) and reaction energies (DG, kcal mol −1 ) given in parentheses for conrotatory electrocyclization of 2-substituted (entries 1-7), (3Z)-methyl substituted (entry 8), (3E)-substituted (entries 9-14), (4Z)-methyl substituted (entry 15) and (4E)substituted (entries 16 for the unsubstituted 1,3,5,7-tetraene.It is worth mentioning that the strength of the electron withdrawing character is directly linked with facilitating the re-arrangement of the pbonds and generation of the new s-bond, to form a ring structure.Overall, a comparative evaluation of the results extracted for C1-to C4-substituted tetraenes with the unsubstituted parent 1,3,5,7-octatetraenes substrate shows a trivial impact on the acceleration of the 8p electrocyclisation reaction. With regard to further verify the suitability of the selected DFT method for this investigation, calculations using the M06-2X/6-31+G(d,p) level of theory are conducted for all the monosubstituted patterns.The computed values (refer to ESI †) express only minor variation by up to 0.5 kcal mol −1 and 0.9 kcal mol −1 for the activation energies and the reaction energies, respectively when compared with M06-2X/6-31+G(d) level of theory.Henceforth, we conclude that the use of a lower basis set which is less computationally expensive is suitable to generate a reliable information for this organic transformation.

8p electrocyclisation of disubstituted 1,3,5,7octatetraene
We have also extensively explored the disubstituted 1,3,5,7octatetraenes using the M06-2X/6-31+G(d) method (Table 4).All possible substitutional patterns of different combinations at different site on the tetraene backbone are computed.Such investigation aims to understand the activation barrier requirements for the 8p electrocyclisation of the conjugated difunctionalised system.
2.3.1 1,X-Disubstituted 1,3,5,7-octatetraenes (X = 2-8).The activation energies for 1,2-disubstituted octatetraene are inuenced by both the electronic and steric effects of the substituent on the molecular backbone and, according to their electronic nature.Inclusion of strong electron donating group (-NH 2 ) at both positions (C1 and C2) on the tetraene skeleton demonstrate a noticeable accelerated effect to the tetraene cyclization, reected by the reduced activation energy to 9.7 kcal mol −1 from 17.0 kcal mol −1 (Table 4).This is related to the capacity of the electron donating character of the substituents on one of the terminal carbon double bonds (C1-C2), in which they enhanced the orbital interaction between the two termini (C1 and C8) to form a ring structure.
In the case of dinitro functionalized system (1-NO 2 -2-NO 2 ), it is not surprising that only minor reduction in the activation barrier by 3.3 kcal mol −1 from 17.0 kcal mol −1 is recorded, in spite of their strong electron withdrawing effect.It arises from the steric congestion between the vicinal nitro groups that slightly destabilizes the transition structure.Greater destabilization leading to a higher energy barrier is presented by 1-CHO-2-OH substituted compound at 21.0 kcal mol −1 due to the unfavourable adjacent side groups clashes.
In 1,3-disubstituted model systems, the strain experience between the motifs becomes less eminent as they no longer feature in vicinal clashing.The reaction energies are seen to slightly be driven by the electronic properties of the substituents due to their minor electronic contribution at C3 position to the terminal double bonds to facilitate the formation of the new s-bond.Hence, a decrease in reaction reactivity is recorded by the 1,3-diamino substituted system (1-NH 2 -3-NH 2 ) to 21.2 kcal mol −1 , in contrast to the corresponding 1,2-disubstituted (9.7 kcal mol −1 ) and also the unsubstituted system (17.0 kcal mol −1 ).Similar trends are also observed for some other 1,3-difunctionalized tetraenes.However, reverse trend is observed when functionalities with high electron withdrawing character are incorporated.An obvious increased in the activation barrier up to 8.5 kcal mol −1 and 9.6 kcal mol −1 from 17.0 kcal mol −1 , are shown by the 1-NO 2 -3-F and 1-NO 2 -3-NO 2 system, respectively.The effect of the functional groups installed at positions C4, C5 and C6 on the tetraene skeleton of the 1,X-disubstituted octatetraene become less signicant during the cyclisation step.The location of these installation inconsequentially delegates effective electronic contribution with the molecular orbital for the newly formed C1-C8 bond.The energy trends observed for 1,7-disubstituted octatetraene are more pronounced than the other 1,X-systems (Table 4).It is the electronic effect of both substituents at these positions that are responsible for the observed accelerated transformation.Direct attachment of these groups on the terminal carbon double bonds promotes the orbital interactions between the forming C1-C8 s-bond.A measurable decrease in the activation energies of 1-R 1 -7-CHO and 1-R 1 -7-NO 2 systems are observed, particularly as a result of the R 2 side group's strong electronic nature and their position on the tetraene skeleton where the electronic contributions are maximized and steric effects are minimized.The most accelerated transformation is observed by the 1-NH 2 -7-NO 2 system, in which the activation energy is 8.0 kcal mol −1 compared to 17.0 kcal mol −1 for unsubstituted system.The synergistic effect between the two substituents with opposite electronic characteristics (-NH 2 having electron donating effect, and -NO 2 having electron withdrawing effect) and the magnitude of these properties effectively promotes the rate of ring closure.Greater orbital interaction with the terminus carbon atoms renders better stabilization to the transition structure, hence favouring the cyclooctatriene formation.Similar reactivity trends of such effects is also observed for 1-CH 3 -7-NO 2 and 1-OH-7-NO 2 substituted systems at 8.6 kcal mol −1 and 10.7 kcal mol −1 , respectively.
In addition, a slightly more enhanced reactivity resulted by this synergistic effect is shown by 1-NH 2 -8-NO 2 functionalized octatetraene.The computed activation free energy for this combination is 7.4 kcal mol −1 .Despite the slight destabilization by the steric crowding between the vicinal functional groups at the transition state, the capacity of the electron donating and electron withdrawing characters, together with the direct attachment to the terminus position leads to a more stabilized and favourable orbital interaction.On the contrary, the highest activation barrier is shown by 1-CN-8-CN system at 21.4 kcal mol −1 , that is governed by the greater steric clashes between the adjacent nitrile groups on the transition structure.The clashing of these substituents outweighs the electronic properties, that lead to the destabilization of the transition state by 4.4 kcal mol −1 from the unsubstituted tetraene.5).As discussed for the monosubstitution Paper RSC Advances systems, location of substituents on the tetraene skeleton are shown to exhibit an essential electronic effect on the rate of the reaction.Therefore, an insignicant variation between the computed results is shown between the studied difunctionalised 2,4-, 2,5-and 2,6 compounds.As opposed to the unsubstituted compound, a measurable decrease in the activation energies of these difunctionalised system are shown to effectively inuenced by the presence of substituents with strong electron withdrawing characters as represented by the bold values in Table 5.
A remarkably lower activation energy of 6.3 kcal mol −1 is computed by the 2-NO 2 -5-NO 2 octatetraene (Table 5).In the case of 2,7-disubstitution pattern, a much more prominent decrease in the reaction energies is observed which is the most accelerated transformation in this investigation.The lowest activation barrier is computed for 2-NO 2 -7-NO 2 at 4.9 kcal mol −1 compared to the unsubstituted parent molecule (17.0 kcal mol −1 ).In general, any substituent attached at these positions (C2 and C7) offer maximized electronic contribution towards the formation of the new terminal s-bond.Thus, a decrease in activation barriers of the difunctionalised octatetraene from 17.0 kcal mol −1 are depicted.In addition, an increase in the reactivity via lower barrier of certain substitutional pattern with the involvement of opposite electronic nature, presumably arises as a result of the synergistic effect between the attached substituents.For example, the 2-CH 3 -7-NO 2 , 2-NH 2 -7-NO 2 and 2-OH-7-NO 2 reveals an apparent decrease in the energy barriers to 6.9 kcal mol −1 , 7.7 kcal mol −1 , Table 5 Activation energies (DG ‡ , kcal mol −1 ) using M06-2X/6-31+G(d) for electrocyclic ring-closure of variously functionalised 2,Y-disubstituted octatetraenes (Y = 3-7)   6) display similar outcomes as exhibited by the 1,2-(Table 4) and 2,3-substrates (Table 5), whereby the mild steric interaction between the adjacent side groups govern the overall reaction.In regards to the position of the attachment on the tetraene backbone, minimal electronic effect is involved in the organic transformation to the corresponding substituted cyclooctatriene molecule.Therefore, an even higher activation barrier is noted for 3-NO 2 -4-NH 2 compound at 22.8 kcal mol −1 than its corresponding 1,2disubstituted (12.4 kcal mol −1 ), and 2,3-disubstituted compound (13.3 kcal mol −1 ).
In the event of 3,5-and 3,6-patterns, only few combinations reveal a noticeable decrease in the reaction rate, as shown when nitro group is/are incorporated (Table 6).This could be explained in terms of the high electron accepting capacity that is strong enough to facilitate the 8p electrocyclic ring closure.On the other hand, the reaction reactivity for the other attachments is shown to be insignicantly affected by the electronic properties, as shown by the trivial variation from the unfunctionalised octatetraene.
The transition structure of the 4,5-dimethyl substitution features relatively shorter interatomic distances between the adjacent functional groups than the 1,2-dimethyl substitution (Fig. 4).In accordance to the position of the attached functional groups, the rate of the electrocyclisation reaction is dominantly controlled by steric effect and considerably unaffected by the electronic characters of the substituents, as demonstrated by the increased in the activation energies from 17.0 kcal mol −1 (unsubstituted molecule).
In the event of the vicinal attachment involving strongest electron donating group (-NH 2 ), similar trend of increasing activation barriers is observed (from 1,2-to 4,5-disubstitution).However, it is as a result of the decreasing electronic effect on the 8p electrocyclic ring closure in regards to their position on the tetraene skeleton.
Moreover, reactivity trend involving -NO 2 substitution (the strongest electron withdrawing group) also shows a general increase in the activation energies.However, an anomaly is identied whereby the 3-NO 2 -4-NO 2 -substituted octatetraene display an opposite effect (decrease in activation energy to 12.1 kcal mol −1 ).Presumably, extra stabilization of the 3,4disubstitution favoring the electrocyclisation arise from the greater resonance structures at the transition states, compared to the 2,3-and 4,5-disubstituted octatetraene transition structures.Further investigation is underway to explore the crucial aspects pertinent to 8p electrocyclization of 1,3,5,7-octatetraene motifs.

Conclusions
Salient conclusions drawn from this computational study towards 8p electrocyclic ring closure of conjugated 1,3,5,7octatetraenes include: (a) An extensive benchmarking studies reveal that M06-2X/6-31+G(d) method is one of the appropriate DFT methods for this investigation.This was supported by the accurate computed value with the measured activation energy (17.0 kcal mol −1 ) of the 8p electrocyclisation of 1,3,5,7-octatetraenes.
(b) Monosubstituted 1,3,5,7-octatetraenes shown to selectively prefer the transition structure with an outward substituent over the inward.The lower steric clashes stabilize the transition structure towards the 8p ring closure process.However, the opposite effects are also observed for some of the mesomerically electron donating groups (-NH 2 and -OH).The SOI between the substituent lone pair and p * C7ÀC8 orbital results in the extra stabilization to the inward transition structure.Therefore, this exceptional case depends strongly on the capacity of the electron donating characters.Amidst the four possible monosubstitution patterns, attachment at C2 position displayed a more pronounced effect on the acceleration of the electrocyclic process due to greater electronic contributions with the terminal carbon double bonds.Therefore, they enhanced the orbital interaction for the generation of new sbond between the C1 and C8 atom.Contrary to the unfunctionalized 1,3,5,7-octatetraenes molecule, a single substitution does not result much impact on the efficiency of the 8p electrocyclic ring closure of studied tetraenes.
(c) In the case of disubstitution systems, attachments on the terminal double bonds renders a general reduction in the activation energies.This is due to the greater electronic interaction from both substituents that stabilizes the transition structure.Particularly, the involvement of strong electron withdrawing group (-NO 2 ) oen proceeds with reduced energy barriers.A remarkable decrease in the activation barrier up to 4.9 kcal mol −1 from 17.0 kcal mol −1 (parent molecule) is observed for 2-NO 2 -7-NO 2 system.Moreover, the synergistic effect between the attached substituents of opposite electronic nature in a disubstituted environment also reveals an apparent accelerated effect on the reactivity of the functionalised tetraene molecule.For example, a decrease in the activation barrier to 7.4 kcal mol −1 is observed in the 1-NH 2 -8-NO 2 -disubstituted octatetraene.
(d) In the event of vicinal substitution (1,2-, 2,3-, 3,4-, and 4,5systems), the eight electron conrotatory electrocyclisation reaction is disfavoured by the destabilization of the transition state, as a result of the sterically congested vicinal functional group arrangements.Greatest decelerated effect is observed by the 4-CH 3 -5-CH 3 substituted octatetraene at 25.4 kcal mol −1 , which is 8.4 kcal mol −1 higher than the unsubstituted compound (17.0 kcal mol −1 ).A much severe steric clashes are observed in the 4,5-pattern over the other vicinal disubstituted patterns because of the less important electronic effect of the substituents, and the greater steric clashes experienced between the adjacent side groups.

Computational methods
Density functional theory (DFT) 22 calculations were performed using Gaussian 16 (revision C.01) 23 and the GaussView (Version 6) 24 was used to generate input geometries and visualize output structures.For the calculations related to the comparative benchmarking study for the electrocyclisation of unsubstituted octatetraene, hybrid three-parameter functional commonly known as B3LYP 25 and metahybrid functional M0-62X 26 levels of theory with different basis sets were used to locate all the stationary points involved.Both B3LYP/6-31G(d) 15 and M0-62X/6-31+G(d,p) 18 levels of theory have already been reported for the 8p electrocyclisation of tetraenes.Geometry optimization for monoand disubstituted octatetraenes were performed using the M0-62X functional with the 6-31+G(d) basis set, followed by the frequency calculations at the same level of theory.All stationary points were characterized as minima or transition structures based on normal vibrational mode analysis.Thermal corrections were computed from unscaled frequencies, assuming a standard state of 298.15K and 1 atm.Representative transition states were also linked to their corresponding minima through intrinsic reaction coordinate (IRC) calculations. 27All-s-transoid conformations of the octatetraenes have been used to calculate the energies of the reactants.For the substituents (-OH and -CHO) exhibiting more than one conformation, lowest energy conformers have been provided here.
Fig. 3 (a) General (generic) atomic orbital interaction in eight-electron conrotatory electrocyclic reactions.Secondary orbital interaction between the lone pair of the inward substituent with the atomic orbital on C7 of the tetraene compound shown in blue.The atomic orbitals interaction of the substituents, (b) -CH 3 , (c) -NH 2 , (d) -OH, and (e) -F, with that of the C7 position are represented in red and green.
2.3.2 2,Y-Disubstituted 1,3,5,7-octatetraenes (Y = 3-7).Similar energy trends are observed for the 2,3-disubstituted with the 1,2-disubstituted system, where they are driven by the clashing of the vicinal moieties.Mild destabilization of the transition structure caused by the steric congestion of the adjacent nitro groups in 2-NO 2 -3-NO 2 compound are recorded at 19.0 kcal mol −1 , that is 2.0 kcal mol −1 higher than the parent molecule (Table